f <- system.file("external/test.grd", package="raster")
r <- raster(f)
names(r)
xyplot(test~y, data=r, alpha=0.5)
## Not run:
# ##Solar irradiation data from CMSAF
# ##Data available from http://www.box.net/shared/rl51y1t9sldxk54ogd44
#
# old <- getwd()
# ##change to your folder...
# setwd('CMSAF')
# listFich <- dir(pattern='2008')
# stackSIS <- stack(listFich)SISmm <- SISmm*24 ##from irradiance (W/m2) to irradiation Wh/m2
# setwd(old)
#
# names(SISmm) <- month.abb
#
# ##Relation between the January & February versus July radiation for four
# ##differents longitude regions.
# xyplot(Jan+Feb~Jul|cut(x, 4), data=SISmm, auto.key=list(space='right'))
# ##Faster with hexbinplot
# hexbinplot(Jan~Jul|cut(x, 6), data=SISmm)
# ## End(Not run)
Run the code above in your browser using DataLab